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Confinement and wall effects are known to affect the kinematics and propulsive characteristics 
of swimming microorganisms. When a solid body is dragged through a viscous fluid at constant 
velocity, the presence of a wall increases fluid drag, and thus the net force required to maintain speed 
has to increase. In contrast, recent optical trapping experiments have revealed that the propulsive 
force generated by human spermatozoa is decreased by the presence of boundaries. Here, we use a 
series of simple models to analytically elucidate the propulsive effects of a solid boundary on passively 
actuated filaments and model flagella. For passive flexible filaments actuated periodically at one 
end, the presence of the wall is shown to increase the propulsive forces generated by the filaments 
in the case of displacement-driven actuation, while it decreases the force in the case of force-driven 
actuation. In the case of active filaments as models for eukaryotic flagella, we demonstrate that 
the manner in which a solid wall affects propulsion cannot be known a priori, but is instead a 
nontrivial function of the flagellum frequency, wavelength, its material characteristics, the manner 
in which the molecular motors self-organize to produce oscillations (prescribed activity model or 
self-organized axonemal beating model) , and the boundary conditions applied experimentally to the 
tethered flagellum. In particular, we show that in some cases, the increase in fluid friction induced 
by the wall can lead to a change in the waveform expressed by the flagella which results in a decrease 
in their propulsive force. 

PACS numbers: 47.63.Gd, 47.15.G-, 87.17.Jj 



I. INTRODUCTION 

In the highly viscous environment inhabited by micro-organisms, locomotion is a difficult task, and one rarely 
achieved in the absence of fellow organisms, boundaries or other obstacles. Because drag forces from the fluid domi- 
nates inertia, swimming becomes a problem for microscopic life qualitatively different from that of larger organisms 
such as fish, and nature has evolved several strategies for solving it. Flagellated organisms such as bacteria and 
spermatozoa utilize the fluid drag anisotropy of slender filaments (flagella) in order to propel themselves through a 
viscous fluid [1]. From a biological standpoint, both prokaryotic and eukaryotic flagella serve the same purpose, to 
propel the organism through the fluid, but from a mechanical standpoint the filaments are quite different. Bacteria 
such as E. coli and B. subtilis actuate passive helical filaments using rotary motors embedded in the cell walls, and 
whose rotation gives rise to propulsion [2, 3]. In contrast, spermatozoa (and more generally, eukaryotic) flagella are 
active filaments. They possess an internal musculature, termed the axoneme, which deforms in a wave-like fashion 
due to the action of molecular motors. These motors generate time- varying and coordinated bending moments along 
the flexible flagellum, giving rise to traveling waves, and propulsion of the cell [1, 4]. In that case, the waveform 
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displayed by the cell is a physical balance between the motor activity, the flagellum elasticity, and the fluid forces. 

Near solid boundaries, the behavior of both types of swimming cells is strongly affected. Since the governing 
equations for inertialess fluid flow are time invariant, the geometry of the system fully defines the hydrodynamics. 
Bacteria swim in circles near a wall, as the chirality of the flagellar rotation induces a hydrodynamic torque on 
the body [5-7]. Boundaries also tend to hydrodynamically attract swimming cells, and as a result the steady-state 
distribution of motile cells strongly peaks near walls [8-13]. Near walls, large arrays of cilia (short flagella) are known 
to synchronize, and display coordinated modes of deformation known as metachronal waves [4, 14-18]. 

One topic of renewed interest concerns the dynamics of spermatozoa in confinement, as relevant to the situation in 
mammalian reproduction [19]. Early theoretical studies considered flagellated cells with waveforms unchanged by the 
presence of walls. In that case, because of increased drag forces, the cells have to increase their work against the fluid 
to maintain their waveforms, and as a result they speed up when near a boundary [10, 20-23]. If instead the cells 
are assumed to work with a fixed power, the presence of a boundary leads in general to a decrease in the swimming 
speed [20, 21]. 

Physically, the speed at which a cell swims is a balance between the propulsive force generated by its flagellum 
and the drag from the surrounding fluid. Near a boundary, both the flow field generated by the flagellum and the 
subsequent propulsion generated are expected to be modified, but in a manner which has not been quantified yet. 

Recently, an experimental investigation was carried out using optical trapping on human spermatozoa to investigate 
the influence of boundaries on force generation. Briefly, spermatozoa cells swimming near and parallel to a cover glass 
(distance about 5 /im) were optically trapped, and then moved to a pre-defined distance of up to 100 /im from the 
glass surface. As the flagellum of the trapped cell continuously beat, the trap power was then gradually attenuated 
until the cell escaped. The magnitude of the propulsive force applied by the cell, equal to the minimum force required 
to hold the cell in place by the optical trap, was found to be decreased by the presence of the glass surface [24]. These 
results indicate that the cells do not maintain their waveforms, as for a cell with a fixed waveform the propulsive 
force would increase near the wall. This experimental result suggests therefore that, for eukaryotic cells, the interplay 
between flagellum elasticity, internal actuation, and hydrodynamics can lead to a non-trivial relationship between the 
environment (here, confinement) and the propulsive force generated by the cells. 

In this paper, we use a series of simple models to examine the propulsive effects of a solid boundary on passively 
actuated filaments and model flagella. Our work aims at capturing the essential physics that describes the geometric 
effects on the body, and builds on previous studies of flagellar locomotion far from external influences, both for passive 
filaments [25-29], and for active flagella [30-34] (see also Ref. [1] and references therein). By considering different 
modeling approaches for the filament actuation, and by quantitatively including the change in viscous friction due to 
the presence of the wall, we predict analytically the change in flagellar waveform, as well as the resulting change in 
propulsive force. We demonstrate that the relationship between the wall- flagella distance and the propulsive force it 
generates is in general non-monotonic. For the case of passive flagella actuated at one end, the presence of the wall 
increases the propulsive forces generated by the filament dynamics for displacement-driven actuation, but decreases 
them in the case of force-driven actuation. In contrast, for active filaments we demonstrate that the manner in 
which a solid wall affects propulsion cannot be known a priori, but is instead a complicated function of the flagella 
frequency, wavelength, their boundary conditions and the manner in which the molecular motors self-organize to 
produce oscillations. 

The paper is organized as follows. We start with a summary of the general class of elastohydrodynamics problems, 
and prescribe our course of action for determining the propulsive force and force gradients in the presence of a wall for 
our models (Sec. II). Following this we consider a passive filament actuated at one end, and the modification of the 
thrust it produces in the presence of the no-slip boundary (Sec. III). We then consider two models for active flagella 
swimming very close to the wall, first the case of a prescribed internal sliding force, and then the more realistic case 
of flagellar beating via self-organization of molecular motors in the axoneme (Sec. IV). In both cases, we determine 
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FIG. 1: (color online) Schematic diagram of a flagellated cell (or passive flexible filament) in the presence of a stationary no-slip 
boundary, (a): 3/4 view of the model system. The planar beating of the nagellum or filament takes place parallel to the plane 
of the wall (x, y), and the cell is located at a distance h from the surface. For simplicity we ignore the head effects in this work; 
this is a reasonable assumption as many experiments to measure the propulsive force would anchor the head in place (e.g. an 
optical trap or micropipet), (b): Side-view of the model system. Note that from this angle the nagellum or filament appears 
only as a straight rod, as it is assumed to beat along the y direction. 

the propulsive force, and how it is modified by the presence of the wall. We finish with a discussion of our results in 
the context of spermatozoa locomotion (Sec. V) 

II. ELASTOHYDRODYNAMICS AND SETUP 
A. Setup 

The physical system that we will investigate is illustrated schematically in Fig. 1. We consider a single flagellated 
cell (or a synthetic device with a nagellum- like filament), for which the nagellum undergoes planar beating a distance 
h from the surface of a solid boundary. The plane of actuation is assumed to remain parallel to the plane of the wall 
(directions x and y in Fig. 1). In the experiment of Ref. [24], when the cells are trapped, their plane of beating is 
parallel to the surface, and thus it is a reasonable approximation to assume that it remains so at different heights. 
We will consider both cases of a passive filament actuated at one end (Sec. Ill) as well as an active filament with 
internally distributed actuation (Sec. IV) 

B. Hydrodynamics 

Since the Stokes equations are time-invariant, the geometry of the system completely defines the fluid dynamics 
at zero Reynolds number. Although the governing equations are themselves linear, nonlinearities in the shapes of 
swimming organisms can make calculating the flow difficult. Furthermore, because the flow is determined by the 
instantaneous shapes of the surfaces immersed in the flow, nonlocal hydrodynamic effects can make the calculations 
impossible to do without numerical analysis. 

Fortunately, for flagellated organisms the slenderness of the swimming surface in question leads to several simplifi- 
cations that can be made to the fluid dynamics. Using the asymptotic limit of a slender filament, the force acted by 
the fluid on the moving filament is approximately given by [35] 



- [C||tt + c±(i-tt)] -v, 



(i) 



where (± and C|| are the force- velocity resistance coefficients derived from resistive force theory (£y < (±) [1], and 
accounting for filament motion perpendicular and parallel to the filament axis respectively; here t is the unit tangent 
vector along the filament, I is the identity tensor, and v is the velocity of the filament as it moves through the 
quiescent fluid. 

In addition to slender ness, for small amplitudes of the flagellar beat the geometric terms simplify considerably. 
Although real flagella beat with a large amplitude, experiments and numerics have shown that corrections to linearized 
dynamical shape equations are sub- leading [36-38]. In the linearized regime, only the normal velocity component of 
the filament is important for calculating the fluid force, as the tangential motion enters only for higher-order curvature 
or bending of the filament. In that case, we can represent the amplitude of the filament perpendicular to the propulsive 
direction as a function y(x,t), and the fluid drag on the filament can be simply expressed in terms of the amplitude 
y as f f i w -C±dy/dte y . 



C. Flexibility and activity 

In order to model propulsion, we need to balance the expression for the fluid force on the filament with the internal 
forces, iintt of the model flagellum. In the case of a passive filament actuated at one end, only the elastic forces 
contribute to this term. For a model of spermatozoa, the internal distribution of molecular motors in the axoneme 
leads to an active bending moment that contributes to force and torque balance. In the inert ialess realm of low 
Reynolds number that is inhabited by the swimming cells that we examine, the total force on each infinitesimal 
element of the flagellum sums to zero, and therefore mechanical equilibrium is written as 

f/Z + Unt = 0. (2) 

Thin passive filaments are well modeled by the elastic beam theory [39] . Their elastic strain energy associated with 
deformation is given by 

L 

K(s) 2 ds, (3) 

where A is the bending rigidity of the filament, L its length, and k(s) its local curvature along the arclength, s. 

For an active filament that is powered by an internal musculature, such as a eukaryotic flagellum, we must consider 
not only similar elastic restoring forces, but also any internal actuation forces. We use in this paper the model of 
Camalet & Julicher [32] to address the mechanics of active filaments, as illustrated in Fig. 2. The active filament 
is assumed to be composed of two inextensible elastic beams which are attached to the basal (or head) region, but 
allowed to bend relative to one another, and acted upon by a distribution of equal and opposite active force f(s). As 
the filaments bend, they induce a distribution of sliding displacements, A(s), given geometrically by 

A(s) = [ an(s')ds', (4) 
Jo 

where a is the fixed distance between the filaments at the base. The work done by the filament against the internal 
forces is then added to the enthalpy functional, and we get 

(5) 

In the case of a passively actuated filament, the internal forces are zero, and thus only the elastic contributions enter 
the equations of motion. Extremizing the energy given by Eq. (5), for a particular form of the active forces, yields 
the total internal force per length on the flagellum, f^, which must then be equal and opposite to the fluid force. 




E 



ds 
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K 2 (s) + fA(s) 
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FIG. 2: (color online) Schematic of a flagellated micro-organism swimming via an active flagellum. The inset shows the model 
of the active filament that we use in this analysis (adapted from Ref. [32]). The two pieces of the filament slide past one another, 
each exerting an equal and opposite active force f(s) that is dependent on the position along the flagellum (the arclength s). 
This produces a sliding displacement A. In the simplest case, that of purely elastic response, and for small amplitude motion, 
this displacement A induces a restoring force that is proportional to the magnitude of f(s). The distance between filaments, 
a, has been exaggerated for viewing purposes. 



D. Propulsion and wall effects 

Given a system, either the passive or active filament, and a set of boundary conditions for the flagellum corresponding 
to a given physical situation, the elastohydrodynamics balance allows us to solve for the flagellar beating pattern and 
thus for the propulsive force. For small-amplitude motion, the propulsive force, defined as the force acted by the 
beating flagellum on the surrounding fluid when it is beating but not swimming, is given by 



(Cii - c±) f 

Jo 



— —dx 

dt dx 



(6) 



where e x is the unit vector along the average position of the beating filament [28] (see Fig. 1). In particu- 
lar, if y(x,t) deforms as a pure traveling wave propagating in the +x direction, y(x,t) = yo(x — ct), we get 

rL, 



F = 



(C-L -C||)c/ (dy /dx) 2 dx 



e x . In that case, the force on the fluid is in the -\-x direction, and therefore 
the force on the filament, and the swimming direction if it was free to swim, is in the — x direction. 

To determine the values of (± and £y in Eq. (6), we look at previous work calculating friction coefficients for slender 
bodies near walls [4]. For a slender body of radius a and length L in an unbounded fluid the drag coefficients can be 
calculated asymptotically in the slender limit a <C and are given by 

(Cl)o ° a M2L/a)+C 1 (C||)o °~ ln(2L/a) + C 2 (7) 

where /i is the fluid viscosity and C\ and C2 are 0(1) constants that depend on the specific geometry of the filament. 
Near a boundary, if h is the distance between the filament and the wall, the situation relevant to the experiments in 
Ref. [24] is that of h < L. In this near-field limit, the resistance coefficients relevant to the planar beating geometry 
considered in Fig. 1 are given by 

The far-field limit, h > X, yields only a small correction to the values of (± and £n in Eq. (7) which would likely be 
too small to be measured experimentally, and thus we will not consider this case here. For a comprehensive account 
of these calculations see Ref. [4] and references therein. Experiments on sedimenting cylinders near boundaries were 



6 



considered in Ref. [40], where it is shown that the difference in drag between a cylinder far from a wall, and one in 
the near-field regime that we consider can easily be as much as 50%. In particular, the result of Eq. (8) implies that 
there is a gradient in the fluid friction if the body changes its distance, h, from the wall. The friction gradient is given 
by 

^ = _^L ( 9) 

which is always negative, reflecting physically that viscous forces are increased by the presence of a solid boundary. 

For a slender rod dragged through a viscous fluid with a fixed velocity, the sign of this force gradient implies that 
moving the rod closer to the wall will increase the force required to maintain its speed. For an actuated filament, 
friction from the surrounding fluid plays a dual role. It first affects the propulsive force, as given by Eq. (6), through 
change in the drag coefficients. In addition, for either boundary- or internally-actuated filaments, the change in (± and 
C|| in Eq. (1) modifies the force balance, Eq. (2), thereby changing the shape of the flagellum, and therefore affecting 
the propulsion in Eq. (6) through a modification of y(x,t). It is this interplay between viscous friction, elasticity, and 
activity that we propose to quantitatively analyze in this paper. As we detail below, it usually results in a non-trivial 
and non-monotonic relationship between wall distance and propulsive force. 



III. PASSIVE FILAMENTS 

We first consider the case of a passive elastic filament driven at one end. In that case, the elastic energy is 

where the linearized regime allows for the approximation of the curvature n by the concavity of the function y(x). 
Calculating the functional derivative of the energy in Eq. (10) leads to the elastic force density which, when balanced 
with the fluid force density, results in the linearized dynamics equation 

dy _ d 4 y 

as obtained in previous studies [25, 28]. We will consider a harmonic driving with frequency lj, and by focusing only 
on post-transient effects, we will assume a similar periodic dynamics for the filament. 

Since the filament is driven only at one end, the boundary conditions at the tail (x = L) end are 

^0OM) = °' (!2) 

= o. (13) 

This guarantees that the ends of the filament are force- and torque- free. At the driving position (x = 0), different 
types of physical actuation could be implemented experimentally, leading to different boundary conditions. If the 
tangent angle is prescribed, with the position of the filament fixed in place, the boundary conditions are 

2/(0, t) = 0, (14) 
dy 

— (0,t) = ecoscjt, (15) 

where e <C 1 is the magnitude of the tangent angle deviation from horizontal and uj the driving frequency (pivoted 
case). Alternatively, for experiments that involve manipulating a passive filament via optical tweezers, the driving 
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position would be made to oscillate harmonically under torqueless conditions, leading to the boundary conditions 

2/(0, i) = y cos tot, (16) 
A0(O,t) = 0, (17) 

where 2/0 is the magnitude of the oscillation in space (tweezed case). Finally, we could also consider the case where 
an oscillating force or torque is applied to the driving end leading to boundary conditions 

2/(0, t) = 0, (18) 
A q^(0, i) = M cos out (Torqued), (19) 

A^(0,t) = F cos out (Forced). (20) 
ox 6 

Here Mo and fo are the magnitudes of the time- varying torque and force, respectively. These boundary conditions 
described above are experimentally realizable, for example using micro-pipets. 

Since we ignore transient effects, the steady solution for the amplitude can be written as y(x,t) = 5ft \y(x)e~ lUJt } . 
Additionally, we can define a natural length scale through the dimensionless "sperm number", Sp = L/£ u where 
l u = (A/Cicj) 1 / 4 ; Sp is the only dimensionless number in this problem, and as such it fully governs the filament 
dynamics. For Sp <C 1, the penetration length is much larger than the length of the filament, and thus oscillations 
will not decay along the length of the flagellum. For a passive filament, this is equivalent to a very rigid rod being 
actuated back and forth. Conversely, for Sp ^> 1, oscillations decay very quickly, indicating a floppy string or very 
viscous fluid. For an active filament, there are additional length scales that must be considered, as will be discussed 
later. 

We nondimensionalize time by 1/gj, the direction x along the filament by L, and the filament amplitude by yo which 
is specified in the tweezed case, and can be related to the forcing parameters in the other three cases. Specifically we 
have yo = eL for pivoted actuation, yo = F0L 3 /A for forced actuation, and yo = M0L 2 /A for the torqued condition. 
By doing so, we obtain the following dimensionless equation for the amplitude y as 

y"" ~ iSp 4 y = 0. (21) 

The dimensionless boundary conditions for the tail are thus given by 

y"(l) = 0, (22) 
y'"(l) = 0, (23) 

while the various possible boundary actuations at the driving end are 

2/(0) = 0, jf'(0) = 1 (Pivoted), 

2/(0) = 1, y"(0) = (Tweezed), 

2/(0) = 0, y"'(0) = 1 (Forced), 

y(0) = 0, y"(0) = 1 (Torqued), 

With a solution to the amplitude equation, Eq. (21), we find the total force exerted on the filament by the fluid. 
Because we consider harmonic actuation, we only look at the time-averaged propulsive force, which is given, in a 
dimensional form, by 
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where I(Sp) is a dimensionless integral defined by 

I(Sp) = Re 



c^ixxVxx VxVx 



(x = 0), (25) 



where an asterisk denotes the complex conjugate, and Re the real part. This expression is a direct consequence 
of the integral for the propulsive force, Eq. (6), being a total derivative when the filament is passive, and thus the 
function depends only on the non-dimensional amplitude y evaluated at the endpoints [25]; the free end x = 1 does 
not contributed because of Eqs. (22)- (23). 

If we then define a scaling function Z(Sp) = I(Sp)/Sp 4 to contain all of the dependence on Sp for the propulsive 
force, we have that the four cases are given by 

(F) = ^uj( ± e 2 L 2 Z(Sp) (Pivoted), (26) 

(F) = \uUvlZ{Sp) (Tweezed), (27) 

(F) = \ " C± ^ L * Z(Sp) (Forced), (28) 

(F) = \^^Z(Sp) (Torqued), (29) 

and each of the scaling functions is unique to the four boundary conditions that solve the amplitude Eq. (21) 

We can now examine how the propulsive force changes as the filament changes its distance h to the wall. Using the 
chain rule, the force gradient is given by 

d /EA 1 2 fd( ±r7 dZ\ 



taw = r y2 °{it z+ ^oh) ^ 

1 2 fd( ± v , A dZ dSpd( ± 



= ry°{it z+ ^ds- P dciit) w 
= r y °^ [ z + 4 Sp ds- P ) > (32) 

where we have used the definition of Sp to take a partial derivative with respect to the resistive coefficient. Since Sp 
is monotonic in the resistive coefficient f_L, and from Eq. 8, we see that £j_ itself is a monotonic function of h, we can 
examine the qualitative behavior of the force gradient by recasting the derivatives with respect to h in terms of the 
dimensionless number Sp. It is important to note at this point that the variation of Sp with the distance between the 
flagellum and the wall is weak, with a scaling Sp ~ [log fo/a] -1 / 4 , and thus even large changes in the distance from 
the wall will produce small changes in Sp. Quantitative details are discussed in Sec. IV A, and for a bull spermatozoa 
cell with Sp « 7 away from a boundary, bringing the cell closer to the wall will lead to a typical increase in Sp of 
about 10%. 

The force gradient is now composed of two terms: the first is due entirely to the changing fluid friction, but the 
second is a more complicated effect that incorporates the elasticity of the filament through the shape change. This 
term reflects the fact that, in the presence of a wall, the elastohydrodynamic penetration length is a function of the 
distance from the wall. Changing that distance changes the essential character of viscous-induced oscillation in the 
filament. Due to the competition between these terms, it is now not necessarily the case that this gradient be negative, 
as we would expect for the case of the rigid rod. From a dimensional standpoint, we can write 



dh 



Z' (33) 



where we have therefore, since d(±/dh is negative (see Eq. 9), 



Z ' = -[Z+\sp^-\. (34) 
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FIG. 3: (color online) Passive filament: Normalized propulsive force, Z, (blue solid line /top curve) and propulsive force 
gradient, Z' (red dashed line/bottom curve) as a function of the dimensionless parameter Sp, for all four boundary conditions. 
Shapes of the filament for different Sp are superimposed, (a): Pivoted case; (b): Tweezed actuation; (c): Torqued condition; 
(d): Forced boundary condition. In the case of displacement-driven actuation, dZ/dh is negative; since £j_ is also a decreasing 
function of h (see Eq. 9), the time-average propulsive force, Eq. (24), is a decreasing function of /i, and the presence of a 
boundary always increase the propulsion of the filament. On the contrary, for force-driven actuation, the sign of dZ/dh is 
positive, indicating that the presence of the boundary decreases the propulsive force. 



To represent the change of the force with the wall distance, we plot in Fig. 3 the dimensionless force, Z (top, blue 
solid line) and force gradient, Z' (bottom, red dashed line) for all four different boundary conditions. Representative 
shapes of the filament over an entire period of oscillation are also shown for various values of Sp. For low values of 
Sp, the filament behaves like a rigid rod, while for larger values the decay length of the actuated filament becomes 
apparent. In both the pivoted and the tweezed cases, there is a maximum propulsive force (blue curve) that occurs 
for Sp > 1: Sp « 4 for the tweezed case, Sp « 2 for the pivoted case. For the torqued and forced cases, the maxima 
occur for Sp < 1, and as Sp approaches zero these conditions leave the linearized regime, and the model breaks down. 

As can be seen in Fig. 3, Z is always positive meaning that the actuation at x = leads to filament dynamics 
pushing against the fluid always in the +x direction. In addition, we see that Z' is always negative for the pivoted and 
tweezed cases (displacement-driven actuation), indicating that the force (or Z) is a decreasing function of the distance 
to the wall, h. The time-average propulsive force, Eq. (24), generated is therefore always increased by the presence 
of a boundary. As such, it is somewhat similar to the increase in the drag on a body driven at a constant velocity 
by the presence of a boundary. Note however that Z' displays a non-monotonic dependence on Sp (and therefore on 
h, since Sp is a monotonically decreasing function of h), which is due entirely to the change in the filament shape 
accompanying the change in height. In contrast, in the torqued and forced cases (force-driven actuation), we see that 
Z' is positive, indicating that the opposite is true, and the propulsive force is now decreased by the presence of the 
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wall. These results are reminiscent of early work showing a similar contrast between two-dimensional swimming with 
fixed kinematics or fixed hydrodynamic power [20-23]. 



IV. ACTIVE FLAGELLA 

We now turn to the case of an active filament as a model for a eukaryotic flagellum, and investigate how the changes 
in hydrodynamic drag induced by the wall and the response of the flagellar amplitude couples to the internal activity. 
As before, we examine force balance, but for an active flagellum we retain the internal forcing term, and the linearized 
dynamics equation becomes [31, 32, 37] 

, dy A d A y df 

Since the total force and moment on the filament must vanish, this imposes boundary conditions on the distal, or 
"tail" , end as 

-40(I,t) + «/(I) = 0, (36) 

-^OM) = 0. (37) 

Regarding the other side of the flagellum (its "head" side), we will examine two different types of boundary 
conditions, either clamped or hinged. For a clamped filament, physically corresponding to a cell that is immobilized 
via micro-pipet, both the head and its tangent angle cannot change, and the boundary conditions are thus 

2/(0, t) = 0, (38) 

g(0,t) = 0. (39) 

In the hinged case, the head remains fixed, but the tangent angle can move torquelessly. This is the situation that 
would take place in the presence of an optical trap, so the boundary conditions become 

2/(0,*) = 0, (40) 



+ a / f( x ) dx = °- ( 41 ) 



o 



If we non-dimensionalize x by L and time by 1 /lj as in the previous section, but additionally scale the magnitude 
of the flagellar beat yo by afoL 3 /A, where fo is the magnitude of the internal force, then the equation for mechanical 
equilibrium becomes 

f'-^=| (42) 
The dimensionless version of the boundary conditions are 

-y"'(l)+/(l) = 0, (43) 
y"(l) = 0, (44) 

for the tail end, and y(0) =0 and one of the following for the head 

y\0) = (Clamped), (45) 

2/"(0) + / f(x)dx = (Hinged). (46) 



o 



In order to consider the response of the flagellum due to changing the distance from a wall, and thus changing the 
friction of the fluid, the mechanism of the axoneme itself must be taken into account, i.e., a model for / must be 
prescribed. We proceed in the next two sections by considering two such models. 
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A. Prescribed activity 

In the first approach, we consider that the active force per unit length takes the form of a prescribed traveling 
wave, i.e. f(x : t) = ?R{foe tkx ~ lu;t } [31], whose frequency (ou) and wavenumber (k) are not modified by the presence of 
a boundary. This enables us to completely specify the filament shape, its propulsive force, and the propulsive force 
gradient, with two dimensionless numbers: Sp and kL. The dimensionless equation of mechanical equilibrium is now 
written as 

y"" - iSp 4 y = ikLe ikLx , (47) 



with boundary conditions 



2/(0) = 0, (48) 
y"'(l) = e* kL , (49) 

y"{i) = o, (50) 



and either of the following 



y(0) = (Clamped), (51) 
y"(0) = - { — ; (Hinged). (52) 

It is then straightforward to calculate the propulsive force on the fluid due to the active filament, which is given by 

{F) = ^ayo^(Sp,kL), (53) 

where the scaling function Y(Sp,kL) is defined as the dimensionless integral 

T(Sp,kL) = Im [ y*^-dx • (54) 
Jo v x 

Here Im denotes the imaginary part. Note that the amplitude of the flagellar beat is proportional to the magnitude 
of the active force, and thus the propulsive force scales quadratically with it. This dependence has been scaled out of 
the propulsive force, and is assumed to remain constant. 

As in the previous section, the force gradient can be calculated using the chain rule and we get 

The main difference between Eq. (56) and Eq. (32) is that the scaling function T is a function of both the elastohy- 
drodynamic length scale parametrized by Sp and an active length scale defined by kL; the interaction of these two 
lengths leads to several non-intuitive results, as shown below. From a dimensional standpoint, we can write 



d , v 1 9 



dCi 



dh 

where, as above, we have defined the dimensionless force gradient as 



T' (57) 



(58) 
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FIG. 4: (color online) Active filament with prescribed activity and clamped boundary conditions: dimensionless force, T, force 
gradient, Y', and gradient of the square of the norm of the force, YY', for various values of (Sp, kL). The corresponding shapes 
are displayed at their representative values on the curves, (a) T (blue, top), Y' (red, middle) and YY' (bottom, black) as a 
function of kL for Sp = 5 (solid) and 9 (dashed); kL runs from 0.1 to 8tt (b) All three again as a function of Sp for kL = 3tt 
(solid) and kL — Stt (dashed); Sp varies from 3 to 10. The axis limits were chosen to cover a wide range of biologically relevant 
filaments and viscosity solutions, keeping in mind that for small Sp the rigid rod limit renders the model inaccurate and is 
irrelevant for biological locomotion. 



Anticipating that the mean force, (F), can change sign, we compute the gradient in the norm squared of the force 
and get 

ClTT', (59) 

so that YY' is the dimensionless gradient in the norm (squared) of the propulsive force. 

In Fig. 4a we display Y (top, blue), Y' (middle, red) and YY' (bottom, black) for Sp = 5 (solid line) and Sp = 9 
(dashed line) as a function of kL in the case of clamped boundary conditions. In Fig. 4b we show the same system, 
only this time several representative values of kL are chosen and the resulting dimensionless force and force gradients 
are plotted as a function of Sp. Similar results are shown in Fig. 5 in the case of hinged boundary conditions. 

First, we observe that although the activity wave is always traveling from the body of the cell to the tip of the active 
filament (uo/k = c > 0), the propulsive force can change sign: Y > means a force on the fluid in the same direction as 
the wave, and therefore (if the cell was free to move) swimming in the direction opposite to the wave. As a difference, 
Y < means the generation of a force on the fluid opposite to the wave propagation, and therefore swimming along 
the direction of the wave propagation. Recall that in the passive case analyzed in the previous section, we always 
had a positive force. We also observe here, in general, a non-monotonic variation of the propulsive force with both 
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FIG. 5: (color online) Active filament with prescribed activity and hinged boundary conditions: Dimensionless force, T, force 
gradient, Y', and gradient of the norm squared of the force, YY', for various values of (Sp, kL). The corresponding shapes 
are displayed at their representative values on the curves, (a) T (blue, top), Y' (red, middle) and YY' (bottom, black) as a 
function of kL for Sp = 5 (solid) and 9 (dashed); kL varies from 0.1 to 8ty (b) All three again as a function of Sp for kL = 3ty 
(solid) and kL = Sty (dashed); here Sp runs from 3 to 10. 



the activity wavelength (through kL) and frequency (through Sp). We further note the importance of the boundary 
conditions as markedly different results are obtained in Fig. 4 and Fig. 5. 

The second important results to note from Figs. 4 and 5 are the variations in the sign of the force gradient. In 
all cases where the product XT' > 0, the presence of the boundary causes the magnitude of the propulsive force to 
decrease. The opposite is true in all cases where XX' < 0. It should be obvious from Figs. 4 and 5 that the sign of 
the force gradient displays a complex dependence on the parameters (Sp, kL), as well as on the type of boundary 
conditions considered. 

To demonstrate the physical significance of our results, we consider as an example the case of a bull spermatozoa. 
For such a cell, we have L « 60 /im, A « 10 -21 Nm 2 , uo « 10 Hz and the "bare" resistive coefficient (C±)oo ~ 10 -3 
Ns/m 2 far from the wall, leading to Spoo ~ 7 [4]. As discussed above, the wall increases fluid drag, and thus it 
increases the value of Sp. Since is possible to change the drag coefficient by as much as 50%, and since Sp ~ C 1/4 , the 
change in the value of Sp can be as high as about 10%. The active length scale kL is more difficult to estimate, because 
the prescribed activity is not immediately obvious through direct observation of the flagellar beat, but reasonable 
estimates give kL between 3tt and 5tt [31, 41]. In the context of the prescribed activity model studied here, kL is 
assumed to not change with distance from the wall. 

We show in Fig. 6 contour plots of the gradient of the norm of the propulsive force, YY 7 , as a function of both 
Sp and kL (left: clamped conditions; right: hinged conditions). The force gradient is positive in the filled contour 
regions; contour lines are 5 x 10 -3 in dimensionless units of force square per unit length. Any given point in the 
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FIG. 6: (color online) Contour plots for the dimensionless propulsive thrust gradient, XT', as a function of Sp and kL. Left: 
clamped boundary conditions (micro-pipette); right: hinged conditions (optical trap). The force gradient is positive in the 
filled contour regions; contour lines are 5 x 10 -3 in dimensionless units of force square per unit length. The arrows indicate the 
increase in the value of Sp taking place as a beating cell is gradually approached to a solid boundary, starting at Sp = 7 and 
for kL = 37T, 57T and 7ty. 

(Sp,kL) plane gives a particular value of the dimensionless force gradient. By bringing the beating flagella closer to 
the wall, the sperm number is progressively increased, and the value of the new force gradient is found by gradually 
moving along horizontal lines in Fig. 6 (which are lines of constant kL). We show in Fig. 6 arrows corresponding to 
this gradual increase starting at Sp = 7 and for kL = 37r, 5tt and 7tt. 

First we observe, again, that the nature of the boundary condition strongly affects the sign of the force gradient, and 
experiments performed using optical trapping should give different results from experiments employing micropipettes. 
Second, we see that three distinct cases are possible depending on the domain crossed by one of the arrows in Fig. 6. 
In the first case, the cell away from the wall is in a region where XX' is negative (white domains in Fig. 6) and remains 
in it during the increase of Sp; in that case, a measurement would lead to a monotonic increase of the propulsive force 
as the flagellum comes closer to the boundary. A second case is the one for which XX' is always positive (for example 
the middle arrow in the left figure, which remains located inside the positive contour plots), in which case the force 
would be measured to be monotonically decreased by the presence of boundaries. Finally, a third situation can arise 
where the arrow crosses the boundary between a region of positive (respectively negative) gradient and a region of 
negative (respectively positive) gradient, leading to a surprising non- monotonic variation of the propulsive force with 
the flagellum- wall distance. 



The prescribed activity model studied in the previous section assumes that the internal activity of the flagellum is 
not modified by the change in the fluid friction near a wall. In a more physically realistic model, the oscillations of 
the flagellum would arise from a self-organized motion of the molecular motors, and thus the internal force generation 
would in turn be a function of the force distribution on the filament (both bending and fluid drag). 

We consider in this section such a model as first introduced by Camalet and Julicher [32]. In this framework, the 
active force is described by a linear response to the filament sliding, / = %A. Spontaneous oscillations of flagellum then 
arise as a self-organized phenomenon for specific values (eigenvalues) of the response function x- If we use Eq. (35) 
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with the linear response relationship / = x^ ~ axdy/dx, then we obtain the following dimensionless eigenvalue 
equation for the filament amplitude 

y""-xy" + iSp*y = o, (60) 

where x — a L 2 x/A. 

The amplitude can be solved as y = Y^t=i Aie qit ; , where the qi are the four solutions to the characteristic equation 
corresponding to the characteristic modes of Eq. (60). For a given value of Sp, we note that there is an infinite 
number of discrete eigenvalues Xn5 as in previous work, we will consider only the eigenvalue with lowest norm as they 
correspond to the lowest degree of motor activity. Non-trivial solutions to Eq. (60) exist only for certain critical pairs 
of parameters (Sp c ,Xc), se "t by the boundary conditions [32, 37] 

4 

J2M ij A j =0, (61) 

1=1 

where the matrix M is defined by the boundary conditions (clamped or hinged, same as in the previous section). 
Since the solution, y, is an eigenfunction for a linear equation, it is known only up to a multiplicative constant, and 
the absolute magnitude for the oscillations can therefore not be obtained. Away from a wall, it has been shown 
theoretically and experimentally that the beating amplitude depends only weakly on nonlinear corrections [36, 37], 
so we will proceed by writing the (arbitrary) amplitude of the filament oscillations as yo, which we assume is the 
magnitude of the first component of the eigenvector (i.e. A\ = yo), and is assumed to remain constant. 
The propulsive force acting from the flagellum on the fluid in this case is given by 

(F) = ^C±vlMSp), (62) 



where 

A(Sp) = Im 

and which can be written more explicitly as 



.Jo 



y*it dx 



A(Sp) = Im 



J o i= i ,=1 



i qjX dx 



(63) 



(64) 



In Eq. (64), the A^s and <^'s are all implicit functions of the response x-> which is in turn a function of Sp. These 
functions are all known, albeit verbose, and thus allow us to explicitly write the force in terms of Sp. 

Near a wall, the fluid friction modifies the fluid resistance coefficient, but the bending modulus A remains 
constant. In order to elucidate the variation of the propulsive force with a change in flagellum- wall distance, we 
thus need to know how both the beat frequency, cj, and the response function, x, vary. Without further biological 
information about the behavior of molecular motors under changing load, we now have to make modeling assumptions. 

The functional dependence of the response function x-> on the oscillation frequency and other materiel parameters 
is, in general, unknown. In order to satisfactorily solve Eq. (60) we need to find the complex eigenvalues that allow for 
non-trivial solutions of the equation to exist. To fully model the active system, the filament response function should 
be derived from a model for the molecular motors, and without such an explicit model, x could very generally be a 
function of on Sp, ATP production and concentration (i.e. activity level), load distribution, structural inhomogeneities, 
and any other parameter (s) that govern activity in the axoneme. The most general solution would thus require a 
detailed model of the molecular motors (see e.g. [32, 36, 37]). 

Since our primary focus in this paper is to explore the hydrodynamic consequences of activity, boundaries, and 
elasticity, we will make below several modeling assumptions in order to examine extreme cases, by essentially specifying 
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FIG. 7: (color online) Dimensionless force (A, blue line) and force gradient (A 7 , red dashed line) for the linear response model 
with clamped (left) and hinged (right) boundary conditions, as a function of Sp. Characteristic flagellar beat patterns are also 
shown at several values of Sp. 



the functional dependence of x- In general both the frequency uo and the response function x will vary, but we are 
going to assume here that one of them remains essentially constant as the wall-flagellum distance is varied. We thus 
assume that one parameter shows a strong variation with ft whereas the other depends only weakly on ft. 

In the first case, the linear response function remains constant, such that the only way for the eigensolution to 
Eq. (60) to have non-trivial values requires the oscillation frequency, cj, to change in such a manner that Sp remains 
constant. Given the definition of Sp, this means that one would observe experimentally a frequency change given 
by uj(h)/uj(h = oo) = Cj_(ft = °°)/ClW> but no waveform variation. In that case, the only change in the propulsive 
force, Eq. (62), would arise from the variation of (± with ft, and therefore one would experimentally measure an 
increase in the magnitude of the force near boundaries. 

The second, more complex, case is one in which the frequency of oscillation of the filament would remain fixed. 
Although there is no experimental evidence that demonstrates that uo does remain constant as the distance between 
the cell and the boundary is changing, it is a reasonable assumption to make in the absence of this information. This 
assumption, however, puts a stringent constraint on the response function x; in order to provide non-trivial solutions 
to Eq. (60) the response function must be only a function of Sp in such a fashion as to match the eigenvalues exactly. 
In this case, varying the distance between the flagellum and the wall, ft, modifies the value of Sp. As a result of 
changes in Sp, both the response function x, and the eigenfunction, will be continuously modified, which will produce 
a non-trivial variation in the propulsive force. As in the previous section, we will write the force gradient formally as 

A'. (65) 

For this scenario, we plot in Fig. 7 the dimensionless propulsive force (A) and force gradient (A') for the clamped and 
hinged boundary conditions respectively, as well as several exemplary beat patterns. With this modeling approach, 
we see that A is always positive. Using the notation of Fig. 1, the flagellum is therefore always pushing the fluid along 
the x direction, and is thus expected to swim in the — x direction. Furthermore, we obtain that the sign the force 
gradient depends on the nature of the boundary conditions. In the clamped case, we get that A' is always negative, 
and therefore the presence of boundaries systematically increases the propulsive force. In contrast, in the hinged case, 
the function A' is seen to be negative for Sp below 7, and positive otherwise. In that case, and similarly to what was 
reported experimentally using optical trapping in Ref. [24] , the measurements would show a decrease in the propulsive 
force near boundaries. 
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V. DISCUSSION 

Biological cells do not swim in a vacuum; the environment itself is what makes swimming possible, and thus we must 
consider characteristics of the surroundings that may modify motility behavior. Since flagellated organisms generate 
propulsion by actuating an elastic filament to do work against their viscous environment, the specific response of that 
environment can be crucial to understanding the overall locomotion characteristics. The particular environmental 
effect that we have studied in this paper is the modification of fluid drag by the presence of a no-slip boundary, and 
the balance between the deformable waveform of the flagellum and the viscous fluid forces that generate propulsive 
thrust. 

For driven filaments, we have shown that fluid drag plays a dual role: not only does it change the propulsion 
generated by a given filament waveform, but it also affects the waveform itself expressed by the beating filaments. For 
passively actuated filaments, the resulting wall effect is a systematic increase of the propulsive force that the beating 
filament imparts on the surrounding fluid in the case of displacement- driven actuation, while a decrease is obtained in 
the case of force-driven actuation. In contrast, for active filaments as models for eukaryotic flagella, the modification to 
the propulsive force depends sensitively on a combination of the flagellar material properties, the boundary condition 
applied to the flagellum, and the manner in which the molecular motors organize to cause oscillation; different values 
of parameters can increase, decrease or even display non-monotonic influence on the cellular propulsive force. 

Using simple scaling arguments, let us finally estimate the expected size of the propulsive force change induced by 
a wall on an active flagella . Let us compare the order of magnitude for the force far away from the boundary to the 
average change in force between near- and far-field. Far from the wall the propulsive force scales as ~ ^?/o(C±)ooY. 
Since our calculation for the force change due to a wall focuses on the near field, we can estimate the average change 
in force as AF ^Lx dF/dh as L, the cell length, gives approximately the spatial range over which the near-field 
matches with the far-field. Using the estimate dF/dh ~ —uy^d^/dfyY' we therefore get 

AF ^ _ L(dUdh)V _ (g6) 
According to Eq. (9), we have d(±/dh ~ —(^/hfi hence the scaling becomes 

AF - ^ f T 'V (67) 



Foe ft/i(Cl)oo V T . 

Since we know that (Cl)oo ~ /i/ln(2L/a) and (± ~ \ij ln(2/i/a), we get the final scaling relationship 



Foo [\n(2h/a)} 2 \ T J h 

From the results obtained above, we observe that Y'/Y ~ ±1. For human spermatozoa, the parameters are L « 40 
/im, a w 0.20 /im. In the experiment of Ref. [24] the near-field measurements get as close as h = 5 \im which leads to 
AF/Fqo rsj ±1. This simple order-of-magnitude calculation shows that the force could be expected to be changed by 
order one by the introduction of a boundary. In the experiment conducted in Ref. [24] the force was measured to be 
reduced by a factor of three, which is consistent with this simple estimate. 
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